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Abstract 



A measurement of the mass of the W boson is presented based on a sample 
of 5982 W ^ ev decays observed in pp collisions at y/s = 1.8 TeV with 
the D0 detector during the 1992-1993 run. From a fit to the transverse 
mass spectrum, combined with measurements of the Z boson mass, the W 
boson mass is measured to be Mw = 80.350 it 0.140 (stat.) it 0.165 (syst.) =b 
0.160 (scale) GeV/c^. 
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The parameters of the gauge sector of the electroweak Standard Model [jT| can be taken 
to be the fine structure constant, the Fermi constant, and the mass of the Z boson, M^, 
all measured to a precision better than 0.01%. Higher order calculations then relate the 
mass of the W boson, Miy, and the weak mixing angle, 6'vi/, to these three parameters, the 
heavy fermion masses, and the Higgs boson mass. Within the Standard Model, a direct 
measurement of M-^ thus constrains the allowed region for the top quark and Higgs masses. 
Alternatively, a precision measurement of the W mass, when combined with other measure- 
ments of sin^6'vi/, provides a test of the Standard Model. The mass of the W boson has 
been measured recently in a number of experiments 0. We present here a new precision 
measurement. 

We have analyzed a sample oiW ^ ev decays resulting from "pp collisions at \fs = 
1.8 TeV. This sample, which corresponds to an exposure of ~ 12.8 pb~^, was collected 
with the D0 detector during the 1992-1993 run at the Fermilab Tevatron collider. Two 
components of the detector are most relevant to this analysis. The central tracking system 
is used to reconstruct charged particle tracks and the interaction vertex. A central and two 
end uranium liquid-argon calorimeters measure the energy flow over a pseudorapidity range 
1^1 < 4.2 i. 

Both W ^ ev and Z e^e decays are used in the analysis. The electrons from these 
decays tend to be isolated and of high transverse momentum, pt- At the trigger level [Q, W 
candidates were required to have an electromagnetic (EM) energy cluster with transverse 
energy = E sin ^ > 20 GeV and to have missing transverse energy 20 GeV. Here 
]^T= — J2i Eti, with the sum extending over all calorimeter cells. Z candidates were required 
to have two EM energy clusters, each with Et > 10 GeV. 

Offline selection criteria were imposed on the EM energy cluster of each electron can- 
didate. The transverse and longitudinal shower proflles of the cluster were required to be 
consistent with those expected for an electron . The energy leakage of the cluster into the 
hadronic compartment of the calorimeter was required to be less than 10%. The isolation 
criterion of the cluster was satisfled by requiring the total energy within a cone of radius 
R = 0.4 1^, centered on the electron direction, but outside the EM core of the shower 
{R = 0.2), to be to be less than 15% of the energy in the EM core. A spatial match of the 
cluster with a central detector track was required. Electrons with cluster position in the re- 
gion between the cryostats (1.2 < |?7| < 1.5) or within 10% of the boundary of a calorimeter 
module in the central region were eliminated from the data sample. 

Having found events with well-identifled, isolated electrons and for W bosons the required 
]^T, kinematic constraints were imposed on the data. The Et^s of each electron in Z events 
and of the electron and neutrino in W events were required to exceed 25 GeV. The neutrino 
Et was equated to the J^t- In addition, the transverse momentum of the W boson, , had 
to be less than 30 GeV/c. These selection criteria yielded 7234 W —>■ eu events with the 
electron in the central calorimeter {\r]\ < 1.2), 366 Z ^ ee events with both electrons in the 
central calorimeter, and 281 Z — >^ ee events with one electron in the central and one in an 
end calorimeter (1.5 < |?7| < 2.5). 

Since the longitudinal component of the neutrino momentum is not measured, the W 
invariant mass cannot be reconstructed. Rather, the mass of the W boson is extracted from 
the distribution in transverse mass, deflned as mf. = 2 \E^\ \EJ^\ (1 — cos(peu), where (p^u is 
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the angle between the electron and neutrino transverse momenta. The electron direction 
is defined using the centroid of the calorimeter cluster and the weighted average of the z 
positions of the hits on the track. The uncertainty in determining this angle leads to an 
uncertainty of 50 MeV/c^ on M^^. Since the absolute energy scale of the EM calorimeter 
is not known with the required precision, the ratio of the measured W and Z masses and 
the world average Z mass were used to determine the W boson mass. The module-to- 
module calibration of the central EM calorimeter was determined to a precision of 0.5%. 
The energy resolution of the central EM calorimeter has been parametrized for this analysis 
as a IE = 0.015 © 0.13//^© 0.4/E, with E in GeV. The sampling term of 0.13//^ was 
measured in a test beam; the constant term of 0.015l[]:[!?^ was determined directly from the 
observed width of the Z resonance. The uncertainty in the EM energy resolution contributes 
a 70 MeV/c^ uncertainty on Mw- 

The EM energy scale of the central calorimeter was determined by comparing the masses 
measured in 7r° — > 77, J/ip ^ e~^e~, and Z — > e~^e~ decays to their known values [§0. 
If the electron energy measured in the calorimeter and the true energy are related by 
Eraeas = « -Etruc + ^5 the mcasurcd and true mass values are, to first order, related by 
"^meas = « mtrue + ^ f- The Variable / depends on the decay topology and is given by 
/ = 2(i?i+i?2) sin2^/2, where 7 is the opening angle between the two decay products and Ei 
and E2 are their measured energies. Figure |l] shows the constraints on the parameters a 
and 5 obtained independently from the 7r°, the J/ip, and the Z data. When combined, these 
three constraints limit a and 6 to the shaded elliptical region. Test beam measurements al- 
low for a small nonlinear term in the energy response, which affects both a and S and alters 
the ratio Mw/Mz largely through the effect on 6, as shown by the dotted line in Fig. |1|. 

Using the measured masses for the observed resonances, the energy scale factor deter- 
mined is a = 0.9514 ±0.0018t{];|][|5^ and the offset is 6 = -0.158 ±0.015l[]:^? GeV, where the 
asymmetric errors are due to possible nonlinearities. The measured offset is consistent with 
that determined from test beam data, and has been confirmed by a detailed Monte Carlo 
study of energy loss in the central detectors. The dependence of the measured ratio of the 
W mass to Z mass on a and 6 may be estimated from 
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It should be noted that the W mass is insensitive to a if 5 = 0. The uncertainty on the 
absolute energy scale results in an uncertainty on Mw of 160 MeV/c^, of which 150 MeV/c^ 
is due to the statistics of the Z data sample. 

The W mass is obtained from an unbinned maximum-likelihood fit of the data to distri- 
butions in rriT, generated as a function of Mw at 100 MeV/c^ intervals by a fast Monte Carlo 
simulation. This Monte Carlo models both the production and decay of the vector bosons 
and the detector response, and relies heavily on experimental data for input. It starts with 
the double differential W production cross section in p-r and rapidity calculated at next to 
leading order using the MRSA parton distribution functions (pdf) |Tl|. The mass of 



the W boson is generated with a relativistic Breit-Wigner line shape, skewed by the mass 
dependence of the parton luminosity. In the simulation, the W boson width has been fixed 
to its measured value, Tw = 2.07 ± 0.06 GeV/c^ 0. The uncertainty on Tw results in an 
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uncertainty of 20 MeV/c^ on Mw- The W decay products are then generated in the W rest 
frame with an angular distribution respecting the polarization of the W. Radiative decays 
are generated at 0{a) according to p!2 |. 

After generation of the kinematics of the event at the four-vector level, the resolutions 
of the detector are incorporated and the energy scales are set. Minimum bias (MB) events 
are used to model the underlying event, mimicking the debris in the event due to spectator 
parton interactions and the pile-up associated with multiple interactions, and including the 
residual energy from previous beam crossings. The relative response of the hadronic and EM 
calorimeters is established by studying Z events. To ensure an equivalent event topology 
between the W and Z events, Z decays in which one electron is in the end calorimeter 
are included in this study. The transverse momentum balance in Z events is given by 
Pt + Pt + pF'^ + ut = —^T, where ut is the underlying event contribution and p^'^'^ is 
the transverse momentum of the recoil to the vector boson. One finds for the average 
\Pt + Pt + -^Tp = IPt'^P + I'^^Tp assuming = k where is the transverse 

momentum of the Z measured from the two electrons. The cross term on the right hand side 
averaged to zero since the underlying event vector is randomly distributed with respect to the 
Z recoil system. Figure ^ shows the distribution of +Pt + -^rP versus \pt\'^- The data 
shows a linear relation between the EM and hadronic energy scale, and yields k = 0.83±0.04. 
The intercept yields the magnitude of the underlying event vector, \ut\ = 4.3 ± 0.3 GeV/c, 
consistent with the value obtained from MB events. The uncertainty on M^y due to the 
uncertainty on the hadronic energy scale is 50 MeV/c^. 

The recoil against the vector boson is modeled by a single jet. The transverse momentum 



of the W is scaled by k and smeared using a resolution of 0.80/\Jpj^(GeV), as obtained from 
our dijet events. The uncertainty on the jet resolution gives a 65 MeV/c^ uncertainty on 
Mw- The event is superimposed onto MB events, which simulates the underlying event. The 
luminosity profile of these MB events is chosen such that the mean number of interactions 
per crossing is the same as for the W data. 

The modeling of the recoil and underlying event are verified and constrained by com- 
paring the Pt of the Z obtained from the two electrons, p^^, to that obtained from the 
rest of the event: —p^'^'^ — Ut- To minimize the contribution from the electron energy res- 
olution, the vector sum of these two quantities is projected along the bisector of the two 
electron directions. Since ut is randomly oriented and has a magnitude ~ , the width of 
the distribution is sensitive to the underlying event contribution while the mean is largely 
unaffected. The rms of the distribution is 4.44 ± 0.18 GeV/c. The sensitivity of the width 
of this distribution to the mean number of MB events that mimic the underlying event is 
determined by varying the number of MB events in the Monte Carlo. The number of MB 
events preferred by the data is 0.98 ± 0.06, consistent with one. The uncertainty on 
from the underlying event model is 60 MeV/c^. 

The energy underlying the electron was obtained from W events by measuring the energy 
deposited in a region of the calorimeter the same size as the electron cluster but rotated 
away from the electron in azimuth. On average, the underlying event adds 205 ± 55 MeV 
to the energy of central electrons and results in an uncertainty on M\y of 35 MeV/c^. 

Detector and reconstruction biases were also modeled in the Monte Carlo simulation. 
In radiative decays, W — > ez/7, the ev mass does not reconstruct to the W mass unless 



6 



the photon is clustered with the electron. Also, radiative decays in which the photon is 
radiated near, but not fully within, the electron cluster can distort the cluster shape causing 
the electron to fail the shower shape cuts. The same considerations apply to radiative Z 
decays and these effects do not cancel completely in the ratio of the masses. Similarly, the 
recoil system may affect the electron identification, especially if it is close to the electron. A 
measure of the event selection biases, due to electron shape and isolation cuts, is obtained 
by studying the projection of the momentum recoiling against the W along the electron 
direction: u\\ = [p^^'^ + ut) ■ Pt- inefficiency in u\\ would cause a kinematic bias for the 
W decay products. The efficiency as a function of u\\ has been determined from the W data 
using the energy in a cone around the electron, which is used to select isolated electrons. 
The efficiency was verified using Z decays. For u\\ values of 20 GeV there is an inefficiency 
of approximately 10%. The error on Mw resulting from the uncertainty in the u\\ efficiency 
is 20 MeV/c^. 

The QCD jet background in the W sample was determined from an independent jet data 
sample to be (1.6 ±0.8)%. Inclusion of this background shifts the mass by +33 MeV/c^. The 
background from Z —>■ e^e~ events in which one electron is not identified has been estimated, 
using ISAJET ||13l, to be (0.43 ± 0.05)%. Its effect on Mw is negligible. The uncertainty in 
the amount of background, and its distribution in transverse mass, gives an uncertainty on 
Mw of 35 MeV/c^. The 1.3% irreducible background due to — *• rz^ — *• ez/i/z/ was included 
in the Monte Carlo simulation. All other sources of background are negligible. 

The distribution in m-r and the Monte Carlo line shape corresponding to the best fit are 
shown in Fig. ^. The mass, extracted from a fit of the 5982 events in the range 60 < m-r < 
90 GeV/c2, is Mw = 80.350 ± 0.140 (stat.) ± 0.165 (syst.) ± 0.160 (scale) GeV/c^. Table | 
lists the uncertainties in the measurement, which used the MRS A pdf. As a consistency 
check, a fit to the distribution in the range 30 < pf. < 45 GeV/c^ was performed to 
extract the W mass. This fit results in a mass 50 MeV/c^ lower than when measured from 
the rriT distribution. The statistical error on this fit is 190 MeV/c^. 

The largest systematic uncertainty, beyond those mentioned above, is due to the modeling 
of the spectrum and the pdf 's. The correlation between the pdf 's and the pj^ distribution 
has been addressed. To study the uncertainty, parametrizations of the CTEQ3M pdf were 
obtained incorporating all available data and with the W charge asymmetry |jl5| data 
points moved coherently by ± one standard deviation, resulting in a maximum allowed range 
of pdf 's. The parameters governing the nonperturbative part of the spectrum |16| were 
varied simultaneously, as constrained by our measured pf spectrum. The resulting variation 
in the spectrum leads to an uncertainty of 65 MeV/c^ on Mw- 

In conclusion, a new measurement of the W mass from a fit to the transverse mass 
spectrum oiW^ev decays has been presented. The W mass is measured to be Mw = 
80.350 ± 0.270 GeV/c^, where all errors have been added in quadrature. 
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Departments of Atomic Energy and Science and Education (India), Colciencias (Colombia), 
CONACyT (Mexico), Ministry of Education and KOSEF (Korea), CONICET and UBACyT 
(Argentina), and the A. P. Sloan Foundation. 
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TABLES 



Uncertainty 


MeV/c2 


Statistical 


140 


Energy scale 


IGO 


Statistical 


150 


Z systematics 


35 


Calorimeter low energy nonlmeanties 


25 


Other systematics 


165 


Electron energy resolution 


70 


Jet energy resolution 


65 


pdi s, p!p spectrum 


65 


Underlying event model 


60 


T~» 1 1 • 11 • 1 T — 111 (T 1 

Relative nadromc and EM energy scale 


50 


Electron angle calibration 


50 


Energy underlying electron 


35 


Backgrounds 


35 


Radiative decays 


20 


■uy efficiency 


20 


Trigger efficiency 


20 


W width 


20 


Fitting error 


5 


Total 


270 



TABLE I. Uncertainties in the W boson mass measurement. 
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FIG. 1. Constraints on slope a and intercept 5 from observed J/ip ^ e~^e~ (dashed-dotted 
line), vr^ 77 (dashed line), and Z — > e'^e~ decays (solid line). The shaded inner contour 
shows the combined result. The dotted line indicates the allowed area when nonlinear terms, as 
constrained by test beam measurements, are included. 
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FIG. 2. Distribution of \p^^ + versus for Z events. 
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FIG. 3. Best fit to the transverse mass distribution. The arrows indicate the fitting range from 
which the W mass is extracted. 
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